Sudden singularities in semiclassical gravity 
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It has been claimed in a recent paper (T) that sudden singularities will survive in semiclassical gravity. This 
issue is here carefully reviewed, pointing out that such conclusion, even if valid under some specific conditions, 
does not stand in other cases. An explicit example is studied in detail to support our statement, reached in these 
other situations, that quantum effects may in fact drastically modify the behavior of sudden singularities. 
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Since the discovery that our universe is expanding in an accelerated way, a number of papers discussing models with 
phantom dark energy have been proposed in order to explain this unexpected behavior (see, for instance, ||2lEl). These 
models have the common feature that future singularities show up, like (big rip and sudden ones, and others) EIU. Those 
could, in some cases, be avoided using semiclassical gravity. Indeed, it has been shown explicitly in a number of papers 
that these singularities are modified when one takes into account quantum effects due to massless, conformally coupled 
fields |l3l[5-9|. In a recent work HI (see also the preprint |10|) it has been claimed however that, quite on the contrary, 
quantum effects produced by a massive conformally coupled field cannot prevent the formation of sudden singularities, 
the reason being that the renormalized energy of the created particles is zero and, thus, that backreaction of the quantum 
effects will not change the evolution of the Universe near the future singularity. In the present paper we review in some 
detail and try to understand the results obtained in 1 1 1. It will be argued that renormalization of the energy density can in 
fact be performed in a different way, which seems reasonable enough and leads as a result to an energy density different 
from zero. Consequently, it can be concluded that there may be a backreaction available to modify the classical evolution 
of our Universe. In support of our conclusion we will here present a specific example with sudden singularities, where a 
sudden singularity is drastically modified when one takes into account quantum effects due to the conformally coupled 
fields, in the renormalization scheme adopted through this paper. 



In this Section we will review, using a standard language, the work done in Ref. lH. The units employed throughout 
the paper aie c — h = Mp = 1, being Mp the reduced Planck mass. In 1 1 1 a model was considered where the sudden 
singularity is described by the following scalar factor |4 |: 



where ts is the time at which the sudden singularity occurs, as is the value of the scalar factor at that time, and n is a 



I. INTRODUCTION 



II. SUMMARY AND DISCUSSION OF PREVIOUS RESULTS 




(1) 



parameter, satisfying 1 < n < 2. When t — > the universe is in the radiation phase a — > 




(fls — 1) , and when 



t ^ ts the universe is in the singular phase a — > a. 
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On the other hand, the modes corresponding to a conformally coupled, massive scalar field, <f>, satisfy the Klein-Gordon 
equation 

X;: + [fc'+™2a'(??)]Xfc = 0, (2) 

where -q is the conformal time, x = ™d x' its derivative with respect to the conformal time. One can see that, in 
the radiation phase (77 — > 0), Eq. (j2ji becomes x'l + ^^Xfe = Oj ™d that, in the singular phase {t — tg), it reduces to 
Xfe + (A:2+m2a2)xfc =0. 

Further, since equation Q is not solvable analytically, in Ref. [1 1 it is approximated by the more simple one 

x'l + ojl{il)xk=Q. (3) 

with 

'^'^^'^^^ I fc2+m2a2, for > Z 

T]c being the time at which the sudden transition occurs. Note that, this choice for the frequency is equivalent to selecting 
as scale factor in Eq. |2]l 

airj) = ( 0' < (5) 

^ 1 Us, for f] > rye, 



which is discontinuous at i] = rjc- 

Eq. ([3]l, with the frequency (|4]), has the solution 



e~*''''/V2A:, for rj < rjc, 

Xkiv) = { (6) 
afce-'"^"/V2S^ + /3fce*"'=VV2S^, for 77 > 

being ak and /3fc Bogoliubov coefficients, and uj^ — -y/fc^ + m?a1. Matching at rjc, one obtains 

Pfc - ^-7=77—^, (7) 

thus, the energy of each mode is 

4 4 

/ \ I i2 

Pp..,fc(m)= ..,1/3,1 =^^^^-^. (8) 
It foUows from this result that the energy density of the particles produced in the process [see Eq. (5. 1 12) of ifTTl l. namely 



Pparim) = „ „ , / dkk'^ppart.kim), (9) 

diverges. What is here actually important to stress is the fact that the energy density could diverge due to the sudden 
transition. Had the transition been not so abrupt — what indeed occurs with the original equation Q — the density of 
produced particles could have remained finite (a clear example of this phenomenon can be found in lil2J . see also 1 13 |) . 

In im the renormalization of Eq. (j9]l is undertaken. To do that the authors use the so-called n-wave procedure 1 14- 
[JFl, apparently in an incomplete way. To wit, this method is used to renormalize the energy density tensor whose first 
component is the total energy density. For the case of a conformally coupled, massive scalar field this is (see, for instance, 
Eq. (9.152) of na): 

1 f°° 1 f°° cD, 

p(m) = Ppar{m) + pvac{m) ^ „ o 4 / dkk'^ ppart,k{m) + „ o . / dkk^ — , (10) 

and corresponds to the energy density of the particles produced in the process plus the energy density of the zero-point 
oscillations of the vacuum. 
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In fact, this procedure is easy to explain: defining p'^\m) = ^Pnk{™n) where pk(jn) = Ppart,k{'m) + the 
renormaUzed energy per mode is given by llT4ll 



prim) = lim 



(11) 



and, applying the method to Pkim) = Ppart,k{'m) + one obtains p^j^^im) — 0, because p^j^\m) = pk{m). 

However, there is another way to renormalize the energy density, which we will follow here. Since this method was 
introduced, at the beginning, for smooth transitions, that is, when the frequency uJk = -y/fc^ + TTT?a'^{ri) is a smooth 
function-what does not happen for the frequency defined in Q-when the transition is abrupt we will not apply directly 
expression ([TT}. In that case, we go back to the general formulae (see, e.g., Eqs. (9.161)-(9.164) of ifTSl or Eqs. (29)-(31) 
of ifTSI ), which have a wider range of applicability. In the case at hand, one just needs to subtract the energy of the 
zero-point oscillations of the vacuum which correspond in the general situation to the term p^"' (m) (see, for instance, 
the paragraph below formula (21) of L14J ). thus obtaining 

pr{m)^PpartA^), (12) 

which, indeed, corresponds to the energy of the produced particles in the fc-mode. As a consequence, the energy density 
of the particles produced remains divergent. But this could be due to the assumption that the transition from the radiation 
phase to the singular one is abrupt. Anyway, what is clear is that, as a result of our procedure, the renormalized energy 
density is no more zero. 



As an interesting remark we note that one can alternatively obtain the same result (12i by employing the adiabatic 
substraction prescription ifTTil . In fact, it is strightforward to check this extreme using Eqs. (3.1), (3.7), (3.11), and (3.14) 

ofLL8|. 

Further to the point, another difference appear with our procedure, stemming from the fact that with the method used in 
Ol there is no vacuum polarization effect. To be specific, let us consider the scalar, massless, conformally coupled case, 
where the well-known renormalized energy density is given by (see, e.g., ll9l [T9l ) 

^ ^ ' 4807r2 ^ 2 ^ 9607r2 ^ ' 

Taking further into account that p'''^" (0) = ^ (p''^" (0) - T"^™ (0)), with 

^"""(0) = T^iH + 12H^H + 7HH + AH^) + —^{H^ + H^H) (14) 
the anomalous trace ||9], one obtains de renormalized pressure 

p-"(0) = -^(i^ + IIH^H + ,HH + Ih^) - ^^H^. (15) 

If one takes the scalar factor to be constant at late times, a{t) = as, then p'"'^" (0) and p''^" (0) will vanish. But, introducing 
the scalar factor as given by ([TJ, one gets 



1 n{n- l){n - 2){as - I) f t 

P K-J-Q^ 

and 



n— 3 



n— 4 



p-(o)~^ "'"-'"":;"'-^' fi-^i . (H) 



14407r2 ajt^s \ t 

which are divergent quantities at t = tg. Consequently, they may drastically change the behavior of the scale factor at 
t = tg. In other words, quantum effects might actually modify the singularity. More precisely, consider the Friedmann, 
ij2 — ^p, and the Raychaudury equation, H = —^{p + p). Near the singularity, one has 



3{as - 1) 



(--1) ~ n(l-f^'' 



2n{n- 1) / t ^ " ^ 



1 - - . (18) 
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From here, proceeding as in ||5], we consider the Friedmann semiclassical equation 



where p is given by ( [T8| , and then look for singular solutions of this semiclassical equation, with the form 

t 



Hit) =H,-C{1- 



(19) 



(20) 



where Hg, C and n' are unknown parameters. Inserting ( |20| into ( [T9| ) and retaining the leading terms, we obtain 

4807r2(as - 1) 



n + 1, Hs = 12Vi0tt 



\ 




(as - ir 

86407r2a2t2 ' 



C 



Hgasin + 1) 



(21) 



Then, since 3 < n' < 4, it turns out that H, and H do not diverge at t = t^, which means that, for this kind of singular 
solutions, the singularity becomes much milder, owing to the quantum corrections. In fact, if one does not take into 
account these quantum corrections, i.e., using Eq. ([T]), one easily sees that H diverges. 

To show even more clearly that quantum corrections may drastically alter, in some cases, the behavior of future sin- 
gularities we will consider, in the next Section, a specific example that can be studied qualitatively using the theory of 
dynamical system. 



III. A MODEL DRIVING SUDDEN SINGULARITIES 



In this Section we study in detail the solutions of the model, inspired in equations (21) and (59) of ||3j, given by 
the equation of state p = — p — (^r^^^ with 7 > and p^, pi positive parameters, which also gives rise to a sudden 
singularity, because the pressure diverges at p = ps (a finite value of the energy density). To simplify the calculations, we 
will set 7 = 3 and ps ~ pi- Then, solving the classical Friedmann and continuity equations 



P 
3' 



-3H(p + p), 



(22) 



one obtains 



Ps 



1 - (48p,)^/' {ts - tf'^ 



p^-ps (48p,)-'/' {ts-t)-""^ 



(23) 



which proves that, in this model, the universe develops a sudden singularity at time t = ts- 

Note however that this model is very different to the one discussed in Section II. Effectively, the model proposed in 
|[T]|4l satisfies the strong energy-condition p > and p + 3p > 0. On the other hand, in this Section we propose a universe 

containing phantom dark energy (p + p = — (^p'^'jpp < 0). 

Considering once again a scalar massless, conformally coupled field, the anomalous trace is given by lH 



T''^"(0) 



1 



28807r- 



:nR + 



1 



57607r: 



(24) 



where R — 6{H + 2H^) is the scalar curvature and G = 2AH^{H'^ + H) the Gauss-Bonnet curvature. In terms of 
the Hubble parameter, Eq. (j24| becomes Eq. ( [T4| . To obtain the energy density, one has to introduce the trace anomaly 
r'''="(0) = p'^™(0) - 3p'''="(0) into the conservation equation p'~'="(0) = -3i7(p'''="(0) + p'''="(0)) to get p'^™(0) = 
— if (4p''^"(0) — T'"'^"(0)). Since the value of T'"'^"(0) is given by Eq. ([14]), this last form of the conservation condition is 
a linear first order differential equation that can be easily integrated and whose solution is given by Eq. ( 13 1. 



Having obtained the renormalized energy density, the semiclassical Friedmann equation is given by Eq. 
terms of the Hubble parameter has the form 



19 1, which in 



i/2 



1 



14407r- 



■{3H^H + HH - -H^ 



1 



28 807r2 



(25) 
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This is a second order differential equation on H and its solutions could differ from the ones obtained form the classical 
Friedmann equation. To show this we write the semiclassical Friedmann equation and the conservation equation as an 
autonomous system, which general solution is a two-parameter family 



H' 


= Y, 


Y' 


_ 1 


^ 2H 


P' 


= 3iJ 



(26) 



or, equivalently. 



1 

2HY 



{H^ -p- m^Y + Y^ 



(27) 



H+t, H = H/H+, Y = Y/Hl, and p 



mi^ being 



Y (p.-p)3' 

where we have introduced the dimensionless variables t 

H+ = ■\/28807r^, while ' denotes derivative with respect to the time t. 

In the same way as it was proven in |8|, in the contracting region {H < 0) the system (26 1 has, for general solution, a 
two-parameter family of future-singular solutions, namely 




0. 



(28) 



In the expanding phase (H > 0), we look for future singular behaviors of the form p = ps — pi(ts — t)'^ with fj, > 0, 
pi > and ts > t. Near t^, from (26 1 and retaining the leading terms only, we get the particular solution 



l-(36p,)V8(i-,_t)i/2- 



and H 



^s{ts - t) 



1 



\{mpsf'\u-tf'^ 



(29) 



what means that both not only H, but also H' , are finite near ps, after the quantum corrections are duely taken into 
account. 

This means that, when one incorporates quantum corrections the singularity becomes softer in the expanding phase 
{H > 0). Note also that, at the singularity time we have H{ts) = 0, consequently, this particular solution does not enter 
the contracting phase. 

On the other hand, the solutions without future singularities are of the same kind as the ones that appear when one 
considers a universe with a big rip modeled by the equation of state p = ojp, with w < — 1, and one takes into account 
quantum corrections (see, for instance, f8|). These solutions are given by a one-parameter family in the contracting 
Friedmann phase plus a particular solution that converges towards the contracting de Sitter one. To prove this, we consider 
the Friedmann phase (the classical Friedmann solution) given, in dimensionless variables, by 



3 



2 {Hs - H)e 



(30) 



where Hs = \Jps- Then, for ^ 1, by linearizing the system (27 1 around the classical Friedmann solution (30l and 
applying the WKB method (see, e.g., p. 53 — 55 of |20]) one gets the following one-parameter family of solutions which 
converges to the classical Friedmann solution, when H 0, 




( 3 

2 {H^~HY 



\ 




for < 0, 



(31) 



being K a free parameter. To prove that there is a particular solution which tends towards the contracting de Sitter phase, 
one must note that (—1,0, 0) is a critical point of (26 1 (in fact, it is the contracting de Sitter solution). Then, linearizing 

3±\ " 



this system around that point, one obtains a matrix with eigenvalues Ai 
proves the statement. 



-3 < and A± 



> 0, what finally 
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As last point, in order to qualitatively study the system, it is convenient to perform the change of variable p 
After this, the semiclassical Friedmann equation becomes 



\H\. 



d 



or, equivalently. 



3£ PP^s 
8 (Ps - -pf 



(32) 



P 



-dpV{p, p) - 2,ep^p , 



(33) 



where V{p, p) = - | [p^il 



, and e = sign(_ff) 



Eq. ( [32| shows that the system is dissipative (resp. anti-dissipative) in the expanding phase H > Q (resp. in the 
contracting phase H < 0) which means that the system loses (resp. gains) energy in the expanding (resp. contracting) 
phase. As a consequence of this, when the universe is in the expanding phase it loses energy and, due to the form of the 
potential V, it rolls down to p = = ^ and enters the contracting phase H < 0, where it gains enough energy to arrive 
to H = — oo in a finite time, exhibiting the behavior described by Eq. (28 i (a detailed account of this process is given in 



10). 

Here one can directly see the main difference between semiclassical and classical cosmology. In semiclassical cosmo- 
logy, due to the quantum corrections, our universe, which nowadays is in the expanding phase, will bounce and will enter 
into the contracting one, where it will develop a singular behavior described by (28). However, in classical cosmology, the 
expanding and contracting phases decouple, that is, the universe cannot bounce. This is due to the fact that the classical 
Friedmann equation prescribes selecting H ~ iv^' ^'^'^ then, once the sign has been chosen, it cannot be changed any 
more. 

In order to further check this qualitative analysis, we have performed an accurate numerical study, as shown in Figs. 1 
and 2. 



4 

t 




Figure I: H{t) and p{t) obtained by integration of 

being Hq — 0.1 and ps = I. Quantum correction drive the universe to the contracting phase, where the universe has the behavior 
described by Eq. l |28| >. 

To finish, the above results lead to the conclusion that, while almost any solution will develop one or more singularities 
in the contracting phase, the corresponding scalar factor and energy density will go down to zero in finite time, as shown 
in the specific example of the two-parameter family of solutions given in Eq. ( [28] l. Note, moreover, that all remaining 
solutions actually form a zero-measure set of the whole set of solutions. That is, they are unstable, in the sense that 
any small perturbation of their initial conditions will give rise to a solution which develops future singularities in the 
contracting phase. This shows that, in the example considered, quantum effects will drastically modify the expansion of 
the Universe, along with other physical consequences associated with the presence of classical sudden singularities. 



26 1 with initial conditions in the Friedmann phase, i.e., with {Hq, | ^-g -go)° ' ^o) 
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Figure 2: Comparison of the derivatives of H{t). The red line corresponds to the derivative of H{t) without quantum corrections 
(we see that it will diverge at finite time). The blue line correspond to the derivative of H{t} after taking into account the quantum 
corrections, as obtained by numerical integration of the system \26\ . It is plain from the plot that quantum corrections modify the future 
singularity. 
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